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I.  INTRODUCTION 


The  uae  of  oross-correlation  for  detecting  and  tracking  the  source  of 
signals  received  at  spatially  separated  receiving  sites  has  received  consider¬ 
able  attention  in  the  research  community  during  recent  years.  It  is  necessary 
to  understand  the  statistical  nature  of  ambiguity  surfaces  and  the  interdepen¬ 
dence  of  the  cells  in  the  surface  under  realistic  operational  conditions  in 
order  to  accurately  evaluate  the  traoklng  accuracies  and  the  detection  per¬ 
formance  achievable  with  cross-correlation.  The  limited  understanding  of  the 
statistical  nature  of  ambiguity  surfaces  is  based  on  the  statistics  of  a 
single  cell  in  the  absence  of  signal  and  noise  power  fluctuations  and  for 
equal  signal  and  processing  bandwidths  (herein  called  matched  containment, 
refs.  1-3)*  The  effects  of  power  level  fluctuations  and  signal  overcontain¬ 
ment,  where  the  processing  bandwidth  is  larger  than  the  signal  bandwidth,  must 
be  quantified  in  order  to  fully  understand  the  statistical  nature  of  ambiguity 
surfaces  under  realistic  operational  conditions.  The  effects  of  signal  over¬ 
containment  have  been  quantified  in  the  absence  of  fluctuations  (refs.  4-5). 

It  is  well  known  that  signal  noise  power  levels  do  fluctuate  during 
reasonable  observation  intervals  (refs.  6-11).  Signal  and  noise  power  level 
fluctuations  and  the  rate  at  which  power  levels  fluctuate  can  adversely  affect 
the  signal-to-noise  ratios  required  to  attain  a  desired  performance  (ref. 

6).  The  effects  of  fluctuations  on  ambiguity  surface  statistics  are 
unknown.  The  study  results  presented  in  this  report  address  the  effects  of 
signal  and  noise  power  fluctuations  on  detection  performance. 

An  ambiguity  surface  is  a  two-dimensional  function,  fD) ,  which  is 
the  sample  magnitude-squared  of  the  normalized  cross-oorrelation  between  the 
observations  received  at  two  spatially  separated  sites  as  a  function  of  the 
relative  time  delay  (t)  and  relative  Doppler  shift  (fD)  between  the  observa¬ 
tions.  The  surface  is  generated  for  a  specific  integration  time  (T)  and  pro¬ 
cessing  bandwidth  (2Wp)  as  shown  in  Figure  1-1.  In  actual  practice,  the 
processing  bandwidth  is  always  larger  than  or  equal  to  the  signal  bandwidth. 
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Since  the  ambiguity  surface  is  usually  computed  digitally,  the  ambiguity  sur¬ 
face  io  quantized  into  cells  of  width  At  seconds  in  the  delay  dimension  and 
Afp  Hz  in  the  Doppler  shift  dimension,  where  At  _<  1/2Wp  and  AfD  <_  1/T.  The 
actual  structure  and  statistics  are  affected  by  power  level  fluctuations  and 
processing  parameters.  The  accuracy  with  which  the  time  delay  and  Doppler 
shift  can  be  estimated  and  the  ability  to  detect  a  signal  is  in  turn  affected 
by  the  statistics  and  structure  of  the  surface. 

There  are  many  types  of  fluctuation  conditions  which  depend  on  the  rate 
at  which  the  fluctuation  processes  can  vary.  The  types  of  fluctuations  are 
bounded  by  very  slow  fluctuation  conditions  and  very  rapid  fluctuation  condi¬ 
tions.  Very  slow  fluctuation  conditions  occur  when  the  signal  and  noise 
powers  are  unknown  but  remain  constant  throughout  the  observation  interval. 

On  the  other  hand,  very  rapid  fluctuation  conditions  occur  when  the  power 
fluctuations  from  sample  to  sample  are  so  large  that  successive  samples  may  be 
considered  independent.  Then  there  is  the  whole  range  of  fluctuation  condi¬ 
tions  between  the  above  two  extremes.  The  signal  model  used  to  study  the 
effects  of  fluctuation  is  presented  in  Chapter  2.  The  detection  performance 
is  described  in  Chapters  3  and  4  for  slow  and  rapid  fluctuations,  respec¬ 
tively.  The  results  are  summarized  in  Chapter  5. 
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2.  SIGXAL  MODEL 


The  signal  model  used  to  analyze  the  effects  of  power  fluctuations  on  the 
sample  magnitude-squared  correlation  coefficient  (MSCC)  is  described.  The 
model  is  used  to  analyze  the  effects  of  slow  and  rapid  power  fluctuations. 

The  fluctuation  model  is  discussed  in  Section  2.1.  The  probability  law 
governing  fluctuations  and  the  resulting  statistics  are  developed  in  Sec¬ 
tion  2.2. 

2.1  Fluctuation  Model 


There  are  many  types  of  fluctuation  conditions  which  depend  on  the  rate 
at  which  the  fluctuation  processes  can  vary.  The  types  of  fluctuations  are 
bounded  by  very  slow  fluctuation  conditions  and  very  rapid  fluctuation  condi¬ 
tions.  Very  slow  fluctuation  conditions  occur  when  the  signal  and  noise 
powers  are  unknown  but  remain  constant  throughout  the  observation  interval. 

On  the  other  hand,  very  rapid  fluctuation  conditions  occur  when  the  power 
fluctuations  from  sample  to  sample  are  so  large  that  successive  samples  may  be 
considered  independent.  Then  there  is  the  whole  range  of  fluctuation  condi¬ 
tions  between  the  above  two  extremes.  The  procedures  used  to  analyze  the 
effects  of  fluctuations  will  depend  on  the  type  of  fluctuation  condition. 

The  fluctuation  model  is  a  generalization  of  the  zero  mean  complex 
Gaussian  signal  model  that  is  used  to  analyze  the  statistics  of  the  sample 
MSCC  (refs.  1-5) .  The  effects  of  fluctuations  can  be  included  by  modeling 
signals  as  a  compared  process  (refs.  6,  7).  In  this  case,  signals  are  mod o led 
as 

s(t)  *  /p(t)  x(t)  (2.1) 

where  x(t)  is  a  zero  mean,  unit-variance,  complex,  stationary  Gaussian  process 
Independent  of  p(t);  p(t)  is  a  non-negative  random  process  called  the  power 
process.  Slow  fluctuation  conditions  exist  when  the  correlation  time  of  p(t) 
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is  much  larger  than  the  correlation  time  of  x(t) ,  while  rapid  fluctuation 
conditions  exist  when  the  correlation  time  of  p(t)  is  much  smaller  than  the 
correlation  time  of  x(t).  Nonfluctuation  conditions  exist  when  p(t)  is  a 
known  constant  in  which  case  s(t)  is  a  Gaussian  process  with  variance  p. 


Let  Z(&)  be  a  two-dimensional  zero  mean  complex  random  column  vector  with 
elements  z-j (<-)  and  z2(Z)  representing  samples  from  channels  1  and  2  at  time 
lTs  for  1  :  Tg  is  the  sampling  interval,  and  T  =  N-jTg  is  the 

observation  interval.  The  cross-covariance  matrix  of  Z(l)  is  defined  as: 

Rz(H,k)  =  E(Z(1)  z'(k)}  (2.2) 


where  E{*}  denotes  statistical  expectation  and  '  is  the  complex  conjugate  of 
the  transpose.  Let  Z(l)  contain  spatially  uncorrelated  noise  under  the  HQ 
hypothesis  and  contain  correlated  signal  plus  spatially  uncorrelated  noise 
under  the  hypothesis.  Then 


Z(t) 


Y(!) 

SsUfi  x(*)  ♦  /n(TT  y(*) 


(2.3) 


where  S(l)  and  N(A)  are  the  independent  two-dimensional  power  vectors  for 
signal  and  noise,  respectively:  X(l)  is  a  two-dimensional,  unit-variance, 
zero-mean,  complex  Gaussian  random  vector  with  P3eJ®3  the  correlation  coeffi¬ 
cient  between  x1 ( l)  and  x2(&);  and  Y(  A)  is  a  two-dimensional  unit-variance, 
zero-mean,  complex  Gaussian  random  vector  with  independent  components. 


The  sample  MSCC  can  be  computed  from  the  sample  auto-correlation 
matrix.  The  two-dimensional  positive  definite  Hermetian  sample  auto¬ 
correlation  matrix  is 


A 


ZU)  z'u) 


(2.4) 
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The  sample  MSCC  Is  the  sample  magnitude- squared  cross-correlation  coefficient 
between  z^d)  and  z2(t)  and  is  given  by 


(2.6) 


The  PDF  of  can  be  derived  from  the  PDF  of  A  by  ( 1 )  performing  the  change  of 
variables  indicated  in  Eq.  (2.6)  and  (2)  integrating  out  the  auxiliary  varia¬ 
bles  a^,  a22,  and  the  phase  angle  of  a12. 


The  cross-covariance  matrix  of  Z(&)  is  sample  independent  when  the  power 
processes  are  stationary.  In  this  case, 


Rz(!,k) 


(  "* 

|  V’n 


(2.7) 


where  Rg  and  RN  are  the  cross-covariance  matrices  of  the  signal  and  noise 
vectors,  respectively.  Combining  Eqs.  (2.2)  and  (2.7),  we  have 


where 


E{S(i)> 


E{N(4)} 


-lii 


Therefore,  according  to  Eqs.  (2.7)  and  (2.8), 


Ro  =  «N  = 


R1  =  Rs  ♦  Rjj 


N1  0 


0  N„ 


Si  +  N1 


1  s 


S2  +  N2 


P1P2  PTe' 


P1P2  PT® 


(2.8c) 


(2.8d) 


(2.9a) 


(2.9b) 


(2.9c) 


where  Pjj  ia  the  average  power  in  channek  k,  P-j  is  the  true  correlation  coeffi¬ 
cient  between  the  channels,  and  6  is  the  phase  of  the  true  correlation  between 
the  channels.  The  true  MSCC  is  defined  as: 


[  Et/S^DSgU)}  ]  2  P2 


(S^XSg+Ng) 


(2.10) 


“*  A  > 


In  the  absence  of  fluctuation,  Sk  =  S^,  Nk  =  Nk,  and 


SNR-  SNR.  fT 

i  2  3 


(SNR1+1)(SNR2+i: 


(2.11) 


If  S-|(l)  and  S2(^)  are  independent, 


v 


eCi^tid  E{^cTn 


(2.12) 


£  [Ei/sJJ)})2  (e{^S2(A)})2  pg 


(s1+n1)(s2+n2) 


(2.13) 


2.2  Fluctuation  Statistics 


The  two-dimensional  power  vectors,  S(A)  and  N(l) ,  for  signal  and  noise 
are  modeled  as  two-dimensional  Gamma  random  vectors.  This  appears  to  be  a 
reasonable  statistical  model  because  it  is  a  generalization  of  the  distribu¬ 
tion  of  the  observed  single-site  fluctuation  processes  (refs.  6-11).  It  will 
be  assumed  that  the  power  processes  in  channel  1  is  independent  of  the  power 
process  in  channel  2  for  both  signal  and  noise. 


The  PDF  of  the  signal  power  process  in  channel  k  is 


«V  * 


(s/3*1  -(hsa/sJ 


na. 

(sk/Msk)  *  rcM^) 


.  S.  >  0 


» sk<0 


(2.13) 


where  MSk  is  the  signal  degrees  of  freedom  in  channel  k  and  Sk  is  the  mean 
signal  power  in  channel  k.  Similarly, 


MN  -1 

(N  )  K  -(MN  N  /  N  ) 

”  mn! - •  '  »ki° 


«V 


(Nk/MNk)  Kr(MN„) 


(2.14) 


,  Nk  <  ° 


where  is  the  noise  degrees  of  freedom  in  channel  k  and  Nk  is  the  mean 
noise  power  in  channel  k.  According  to  Eq.  (A. 3)  of  Appendix  A,  the  mean  sig¬ 
nal  power  and  variance  of  the  signal  power  is 


\  * 

\  *  <v2/m\ 


(2.15) 


and  for  the  noise,  we  have 


(2.16) 


<V2/MHk 


Therefore,  the  mean  powers  are  unbiased  and  the  variances  vanish  with  increas¬ 
ing  degrees  of  freedom. 


Finally,  the  true  MSCC  is  biased  in  fluctuation  conditions.  A  larger  SNR 
is  needed  under  fluctuation  conditions  to  achieve  the  same  as  for  no  fluc¬ 
tuations.  Substitute  Eqs.  (2.13),  (2.14),  and  (A. 2)  into  Eq.  (2.10).  Then, 


r(MS1  +  1/2)  T(MS2  +  1/2)  SNR^NRg  g 

/ms ^2  f(MS1 )  T(MS2)  (SNR1+1)(SNR2+1)  ps 


(2.17) 
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where  SNRk  =  S^/N^  By  comparing  Eqs.  (2.11)  and  (2.17),  the  bias  factor  is 


BIAS 


r(Ms1  ♦  1/2)  T(MS2  +  1/2) 

/hs1ms2  hms,)  t(ms2) 


\2 


(2.18) 


Eq.  (2.18)  is  plotted  in  Figure  2.1  for  MS^  =  MS2  =  MS.  It  is  compared  to  p| 
for  no  bias.  It  can  be  seen  that  the  BIAS  can  require  significant  increases 
in  SNR. 


3.  SLOW  FLUCTUATION 


Slow  fluctuation  occurs  when  the  correlation  tines  of  the  signal  and 
noise  power  processes  are  auch  larger  than  the  observation  interval.  In  this 
case,  the  power  processes  beoones  an  unknown  constant  but  unknown.  Therefore, 
the  slow  fluctuation  oase  becoaes  the  case  of  constant  but  unknown  signal  and 
noise  power  levels.  The  cuaulative  distribution  function  (CDF)  of  the  sample 
magnitude- squared  correlation  coefficient  (MSCC)  is  derived  in  Section  3.1. 

The  detection  performance  of  the  sample  MSCC  is  presented  in  Section  3.2.  The 
results  are  summarized,  and  the  implications  discussed  in  Section  3.3. 

3.1  Cumulative  Distribution  Function 


The  cumulative  distribution  function  (CDF)  of  the  sample  MSCC  for  known 

SMA's  ia 


P^ptlp1  »P2'Ps,*T^ 


2  2  MT 

pt(i-o,o2o3) 

N--2 


ij,  <'"pt>Vi(vi-ii';pip2ps) 


(3.1) 


where  p§  la  the  threshold,  ps  is  the  correlation  coefficient  of  the  signal 
components ,  Nj  is  the  degrees  of  freedom,  2^1  (• »•?•»•)  is  the  hypergeometric 
function,  and 


SMRk 

°k  *  SNRjg+1  (3*2) 

ia  the  ratio  of  the  SNR  in  channel  k  (ref.  1).  For  slow  fluctuations,  the 
SNR's  are  unknown  constants.  Therefore,  the  CDF  of  the  sample  MSCC  for  slow 
fluctuations  becomes 


P(ptlpa’V  8  f  I  F(ptlpi»p2»ps,NT)f(p1,p2)  dpidp2 


(3.3) 


where  f(p^,P2)  is  the  joint  probability  density  function  (PDF)  of  the  p^'s 


It  is  reasonable  to  assume  that  and  P2  are  independent  because  the 
acoustic  propagation  conditions  to  the  two  receivers  are  different.  There¬ 
fore,  according  to  the  fluctuation  model  discussed  in  Chapter  2, 


r<V  • 


rCMSjj  .  my 

r(Msk)  ro«k) 


<V 


MSp-1  MN^-1 

MS,.  Pk  (1-C^) 


(l.(Vl)Pk) 


'  “iOlci  ' 


,  otherwise 


where 


(3.4a) 


(3.4b) 


MS^  is  the  signal  fluctuation  degrees  of  freedom  in  channel  k,  MN^  is  the 
noise  fluctuation  degrees  of  freedom  in  channel  k,  is  the  mean  signal  power 
in  channel  k,  and  Nk  is  the  mean  noise  power  in  channel  k.  The  CDF  of  the 
sample  MSCC  is  obtainable  by  substituting  Eq.  (3.4)  into  Eq.  (3*3)  and  evalu¬ 
ating  the  integral. 


Eq.  ( 3 - 1 )  becomes,  upon  expanding  the  hypergeometric  function, 


N-2 


F(pt  lprp2’cVMT) 


4=0 


p=0 


(p!)‘ 


t(p | P^ , P2) 


(3.5a) 


where  (x)n  =  r(n+n)/r(x)  is  Pochammer's  symbol, 

-  N 

^(plp^Pg)  s  (PiP2)P(1"p1p2ps)  (3.5b) 

Define 

t(p)  =  f  f  t(p|p1,p2)f(p1)f(p2)dp1dp2  .  (3.6) 
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Under  the  H0  hypothesis,  P{  *  0,  and  the  CDF  becomes 

«»9  v.  * 


.-2 


*  P- 


q=0 


g 


2  V1 

t  -  o-p|)  T 


(3.10) 


Bq.  (3*10)  is  the  same  as  the  equation  for  the  CDP  of  the  sample  MSCC  under  HQ 
for  known  noise  powers,  Bq.  (4.8)  of  reference  3.  This  is  not  surprising 
because  the  powers  normalize  out  of  the  sample  MSCC  whenever  the  noise  powers 
are  constant,  event  if  they  are  unknown  Eqs.  (2.4-2.6). 


I 


3.2 


Detection  Performance 


The  detection  performance  is  quantified  by  the  probability  of  false  alarm 
(PFA)  and  the  probability  of  detection  (Pp).  These  are  defined  as 

pD  *  1  -  ^Ityi  (3.1D 


and 


Pi 


*  V1 

i  -  F(p||wT>0  .  C1-.P*)1 


(3.12) 


These  equations  are  evaluated  for  equal  channel  conditions  where 


MS 

3  MS,  3 

ms2 

(3.13a) 

MN 

3  MN,  3 

mn2 

(3.13b) 

SNR  s  SNR.  s  SNR. 


(1.13c) 


The  performance  is  quantified  by  (1)  solving  Eq.  (3.12)  for  the  threshold,  p|, 
for  a  specified  PFA,  and  (2)  numerically  solving  Eq.  (3.11)  Tor  the  SNR 
required  for  the  specified  PD,  NT,  and  p|. 

The  performance  is  plotted  as  a  function  of  Nj  in  Figure  3>1>  It  is 
immediately  apparent  that  (1)  uncertainty  in  signal  and  noise  powers  can 
require  large  increases  in  SNR  to  achieve  the  same  performance  in  the  absence 
of  fluctuation,  and  (2)  fluctuation  effects  decrease  as  the  uncertainty 
decreases  (i.e.,  MS  and  MN  increase). 

The  performance  is  plotted  as  a  function  of  PD  for  various  NT's  in 
Figures  3. 2-3. 4.  It  is  apparent  that  fluctuation  effects  decrease  as  MS  and 
MN  increase.  SNR  is  more  sensitive  to  signal  power  fluctuations  than  to  noise 
power  fluctuations  because  the  required  SNR  is  larger  for  NS  =  2,  M  *  10,  than 
for  MS  s  10,  MN  =  2.  This  effect  is  larger  for  larger  Pp's  (>  0.4).  The  SNR 
sensitivity  to  signal  fluctuations  follows  from  the  fact  that  the  PpA 
threshold  is  independent  of  noise  power  fluctuations. 

3.3  Discussion 

The  cumulative  density  function  (CDF)  of  the  sample  MSCC  was  derived  for 
slow  fluctuations.  The  CDF  of  the  sample  MSCC  is  independent  of  noise  power 
fluctuations  under  the  HQ  hypothesis. 

It  is  observed  that  the  SNR  required  to  achieve  the  desired  operating 
point  decreased  as  the  fluctuations  decreased.  The  SNR  is  more  sensitive  to 
signal  power  fluctuations  than  to  noise  power  fluctuations  because  the  PpA 
threshold  is  independent  of  noise  power  fluctuations.  Slow  fluctuations  can 
require  a  4-6  dB  increased  SNR  over  the  SNR  required  to  achieve  comparable 
performance  in  the  absence  of  fluctuation. 
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4.  RAPID  FLUCTUATION 


Rapid  fluctuation  occurs  when  the  power  level  fluctuations  from  sample  to 
same  are  so  large  that  successive  samples  may  be  considered  independent.  The 
cumulative  distribution  function  (CDF)  of  the  sample  magnitude-squared  corre¬ 
lation  coefficient  (MSCC)  is  derived  in  Section  4.1,  and  the  detection  per¬ 
formance  of  the  sample  MSCC  is  presented  in  Section  4.2.  The  results  are 
summarized  and  the  implications  discussed  in  Section  4.3. 


4.1  CDF  of  the  Sample  MSCC 


The  CDF  of  the  sample  MSCC  for  rapid  fluctuation  is  difficult  to  derive 
because  the  probability  density  function  (PDF)  of  the  observation,  Zj  defined 
in  Eq.  (2.3),  is  unknown  for  Gamma  distributed  fluctuations.  This  problem  is 
overcome  by  using  an  Edgeworth  series  approximation  to  the  PDF  of  (Appen¬ 
dix  B).  The  only  way  to  incorporte  the  Edgeworth  series  into  the  derivatives 
of  the  CDF  of  the  sample  MSCC  is  to  use  an  Edgeworth  series  approximation  to 
the  PDF  of  the  sample  auto-covariance  matrix  (Chapter  2).  The  Edgeworth 
series  approximation  to  the  CDF  of  the  sample  MSCC  for  rapid  fluctuation  is 
derived  in  Appendix  C.  The  CDF  of  the  sample  MSCC  with  signal  present  is 


G(PtlNT)l  *  F(P?|HT>1  *  V'^Vl 


(4.1) 


where 


Fd^lVl 


N„ 


Nt-2 


(1-p*)  T  £.  (1-p^Vl(MT’k+1;1;PtPT) 


(4.2) 


is  the  CDF  of  P2  for  no  fluctuation  (ref.  1);  p|  e  (o,1)  is  the  threshold;  and 
F(P^|N<r)1  is  the  Edgeworth  series  correction  factor. 


Since  G(pf  |NT)  and  F(P2|nt)  are  CDFs,  F(  1 1 NT)  -|  =  0  because  G(l|NT)i 
F(l|NT),  =  1.  The  correction  factor  is 


N+>  CM  *> 


|Nt)q  =  PR1F(p;|0fNT-2,NT,NT,1) 

+  PB2F(p^lO,NT-2,HT,HT+1,1) 

+  PR3F(p2|0,NT-2,NT,NT+2,1) 

+  PR4F(p2|0,NT-2>NT+1,NT+1,1) 

+  pr5f(p^|o,nt-i,nt+i,nt+i,i) 

♦  PR6F(p2|0,NT-1,HT+1fNT+2,1) 

♦  PR7F(p2|0,NT,HT+2lNT+2,1) 

+  pr8f(p^i»nt-2»nt+1,V1,1) 

+  PRgF(p2|l,NT-2,NT+1fNT+1,2) 

+  pR10F(ptl1,V2,V1,NT+2,2) 

+  PR , ! F ( p2 | 1 , NT-2 , Nt+2 , Nt+2 , 2 ) 

+  PR12F(p2|2,NT-2,NT+2,NT+2,3) 

|o,e,0,4>,y)  s 

p2(oh-i)  a 

Tgnr y^(Jc+1) Jl-P?)k,F  (0,<>lO4.k+1;Y,OM.a+2;p2p2) 


is  Pochammer's  symbol,  where  the  PR's  are 


(Vi  >  3  V2 

PR1  =  CM7  N„+1  <1“pt) 


(N--1 r  5  M_-1  N_(N_-1) 

PR2  =  CN5 Hhr  (1-4>  '  4CN7  nI+1  (1-p 


PR 


N 

3  =  (CN1(1-pf)  T  -  CN5(1-pJ)  T  ♦  2CN7(1-p^) 


_  H-1 
,2>  T 


,-.2, 


PR4  =  (cn4(i-p‘)  -  CN5(1-p‘) 

,  V2  MNT-1) 

-*■  2CN?(1-P2)  ) 


Nt+1 


id  2  V 

PR5  =  ("  CN5  *  2CV  Vi0'# 

2  V 

♦  2CN7Nt(1-P‘) 


PR6  =  (-  2CN,  +  CN3  -  2CN4)  Nt(1-PP)Nt 


N„ 


♦  (CN5(3-P2>  -  4CN6  -  4CN?)  N^l-p2)  T 


PRy  s  (CH1  +  2CN2  -  CN3  ♦  CN4(1+P7) 


-  CN5  +  2CN6  +  CN7)(Nt+1)(1-pJ) 


2  "T 


J.  ,  N  -1 

PR8  =  (CN^I-p*)  T  -  2CN6(1-p‘)  T 

2  2  W” 


*  *evf(H$>  )  -V 


,+1 


,  2  V1 

PR9  =  (-  2CN6(NT-1)(1-pJ)  T 

,  p  N--2  N-(N--1 ) 

+  4CN_N_p,(1-p,)  1  )  ---  , — 

7  TKT  WT  }  Nt+1 


N  N  - 

PR10  =  ("  2CH3(1-p^)  T  +  (2CN5p^  +  CM6)(1-4>  T 

3  3  V2 

-  8cn7p‘(i-p‘)  t  )  ht(ht-d 


PR 
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=  _  (-  (l»CN2  -  C«3  ♦  2CNJ4pJ)(1-pJ) 


2vNT 


(2CN5p£  ♦  2CN6(1-2pJ)  -  4CN7p^)(1-pP) 


nt(nt+i ) 


PR12  s  (CNgd-P^)"1  -  CV^I-p^"7"1 


H  —2 

CH7p^(1-p2)  T  )Nt(Ht+1)(Nt-1) 


where  the  CN's  (Appendix  D)  are 


SNR?  ♦  MS./MN,  SNR?  ♦  NS./MN. 

«1  s  — 3 - —r  +  — 1 - Hr 

1  MS  (SNR  +1) 2  MS  (SNR  +1)2 
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CM, 


BIAS 


(l  -  BIAS1/2) 


SNR 


SNR. 


CN. 


2  Oil 

2pt^ms^Tsnr^TT  +  ms2(snr2+d^ 


CN„ 


PT 


BIAS 


(1  -  BIAS) 


CN. 


CN3  -  2CN1  -  2CN4 


CN, 


CN, 

2CN„  -  -r*  + 


cv? 


CN_ 


CN1  ♦  2CN2  -  CN3  +  CN4(1+p‘) 


and  where 


2  SNR^SNR2  2 

PT  =  BIAS  (SNR1+1)(SNR0+1 )  Pa 


BIAS 


T(HS1  ♦  1/2)  T(MS2  +  1/2) 


/MS HMS^  r(MS2) 


SNRU  a  sk/Nk  ,  and 


(4.6b) 


(4.6c) 


(4.6d) 

(4.6e) 

(4.6f) 

(4.6g) 


(4.7a) 


(4.7b) 


(4.7c) 


P9  is  the  correlation  coefficient  between  the  signal  components. 


Under  the  Ha  hypothesis ,  p|  a  SNRk  =  0.  Then  the  CDF  of  the  sample  MSCC 

becomes 


0(»'|Ht)o  .  F(p*|»  )  ♦  »J(p2|N  ) 


A 


(4.8a) 


3  r  4ht  ♦  2nt-6 

K'Vo  *  c"i|(vi)(NT-i)F(‘VllT), 


(SN.  *  10)  - 

_  4  *  /  l 


VTi-  FW' V'\ 


4P(p||NI«2)< 


,  and 


CN  -  1 -  +  3 

1  '  MN1  MH2 


4.2  Detection  Performance 


(4.8c) 


(4.8d) 


The  detection  performance  is  quantified  by  the  probability  of  false  alarm 
(PpA)  and  the  probability  of  detection  (PD).  These  are  defined  as: 


fD  *  ’  -  °«>tlVi 


(4.9) 


pFA  =  i-wllVo 

Where  Q(p||Nt)1  and  G(p||Nt)q  are  defined  in  Eqs.  (4.1) 
tions  are  evaluated  for  equal  channel  conditions  where 


(4.10) 


-  (4.8).  These  equa- 


MS  =  MSV  =  MS2 


(4.11a) 


Ml  :  MX.  =  MN0 


(4.11b) 


The  equation  expressing  the  relationship  between  the  PpA,  p|,  and  NT  for 
no  fluctuation  is  well  known  (ref.  1).  It  is 


2  V1 

PFA  =  (1  Pt) 


(4.13) 


By  comparing  Eq.  (4.13)  to  Eq.  (4.10),  it  is  apparent  that  rapid  fluctuation 
affects  the  PpA  threshold,  p^.  One  of  the  attractions  of  using  the  sample 
MSCC  for  detection  is  that  is  independent  of  the  noise  properties  in  the 
absence  of  fluctuations.  However,  this  property  does  not  hold  when  rapid 
fluctuations  are  present  in  the  noise.  The  rapid  fluctuation  thresholds  for  a 
specified  PFA  and  N^  are  plotted  in  Figure  4.1  for  various  fluctuation 
parameters.  The  p|  are  computed  by  numerically  solving  Eq.  (4.10)  for  speci¬ 
fied  PpA,  Nj,  and  MN.  It  is  seen  that  rapid  fluctuation  has  the  largest 
effeot  on  q\  for  6  <  NT  <  500.  It  is  also  apparent  that  the  influence  of 
rapid  fluctuation  decreases  with  increasing  fluctuation  degrees  of  freedom 
(MN)  because  the  variance  of  the  fluctuation  process  decreases  as  MN 
increases.  Eg.  (2.16). 


The  performance  is  plotted  as  a  function  of  NT  in  Figure  4.2.  It  is 
apparent  that  (1)  rapid  fluctuation  can  require  large  increases  in  SNR  with 
respect  to  the  SNR  in  the  absence  of  fluctuation,  (2)  fluctuation  effects 
decrease  as  MS  and  MN  increase,  and  (3)  fluctuation  effects  decrease  as  NT 
increases.  This  means  that  the  performance  becomes  somewhat  insensitive  to 
rapid  fluctuation  for  large  NT. 

The  performance  is  plotted  as  a  function  of  PD  for  various  NT's  in  Fig¬ 
ures  4.3  through  4.5.  It  is  apparent  tht  fluctuation  effects  decrease  as  MS 
and  MN  increase.  SNR  is  more  sensitive  to  noise  power  fluctuations  than  to 
signal  power  fluctuations  because  the  required  SNR  is  larger  for  MS  =  10, 

MN  s  2,  than  for  MS  =  2 ,  MN  =  10.  This  sensitivity  follows  from  the  fact  that 
the  PpA  threshold  is  dependent  on  the  noise  power  fluctuations.  It  is  also 
seen  that  the  effects  of  rapid  fluctuation  can  be  reduced  by  increasing  N<p. 
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4.3  Discussion 
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The  cumulative  density  function  (CDF)  of  the  sample  MSCC  wns  derived  for 
rapid  fluctuation.  Due  to  the  mathematical  difficulties  inherent  in  the  deri¬ 
vation,  the  CDF  was  approximated  with  an  Edgeworth  series. 

The  CDF  of  the  sample  MSCC  under  HQ  is  dependent  on  the  noise  power 
fluctuation  parameters.  Consequently,  the  sample  MSCC  loses  some  of  its 
attractiveness  as  a  detector  because  the  PpA  threshold  is  dependent  on  the 
fluctuation  parameters.  This  is  in  contrast  to  the  CDF  of  the  sample  MSCC  in 
the  absence  of  fluctuations  where  the  CDF  is  only  dependent  on  the  noise 
degrees  of  freedom. 

It  is  observed  that  the  SNR  required  to  achieve  the  desired  operating 
point  decreased  as  the  fluctuation  decreased  (i.e.,  MS  and  MN  increased).  The 
SNR  is  more  sensitive  to  noise  power  fluctuations  than  to  signal  power  fluctu¬ 
ations  because  the  PpA  threshold  is  affected  by  the  noise  power  fluctuation. 
Rapid  fluctuations  can  require  a  3-4  dB  increase  in  SNR  over  the  SNR  required 
to  achieve  the  comparable  performance  in  the  absence  of  fluctuations. 
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5.  CONCLUSIONS 


A  detailed  analysis  of  the  detection  performance  of  the  sample  magnitude- 
squared  correlation  coefficient  (MSCC)  in  the  presence  of  fluctuations  has 
been  presented.  Fluctuations  can  be  characterized  as  slow  or  rapid,  or  as 
anything  in  between.  The  fluctuation  process  samples  are  completely  corre¬ 
lated  over  the  observation  interval  for  slow  fluctuation,  while  the  fluctua¬ 
tion  process  samples  are  completely  uncorrelated  over  the  observation  interval 
for  rapid  fluctuations.  These  two  bounds  on  fluctuations  can  be  studied 
analytically.  Simulation  is  required  to  study  fluctuation  processes  with  cor¬ 
relation  times  that  lie  between  the  bounds. 

It  is  concluded  that  fluctuations  require  a  3-4  dB  increase  in  SNR  for 
rapid  fluctuations  and  a  4-6  dB  increase  in  SNR  for  slow  fluctuations  over  the 
SNR  required  to  achieve  comparable  peformance  in  the  absence  of  fluctua¬ 
tions.  In  all  fluctuation  cases,  the  required  SNR  decreases  as  the  fluctua¬ 
tion  becomes  "less"  random  (i.e.,  the  variance  decreases).  However,  the 
effects  of  slow  fluctuation  are  basically  independent  of  the  signal  time- 
bandwidth  product  (Nt),  while  the  effects  of  rapid  fluctuation  can  be  reduced 
by  increasing  N^. 

The  PFA  threshold  (p|)  is  independent  of  the  noise  fluctuation  process 
for  slow  fluctuations,  but  it  is  dependent  on  the  fluctuation  process  for 
rapid  fluctuation.  It  can  be  concluded  that  is  dependent  on  the  fluctua¬ 
tion  process  for  all  correlation  times  except  for  slow  fluctuation.  The 

p 

dependency  of  on  the  noise  fluctuation  process  decreases  as  the  correlation 
time  of  the  fluctuation  process  increases. 

The  SNR  for  rapid  fluctuations  is  more  sensitive  to  noise  power  fluctua¬ 
tions  than  to  signal  power  fluctuations  for  all  fluctuation  processes  with 
correlation  times  less  than  the  observation  interval.  This  is  caused  by  the 
fact  that  (1)  is  dependent  on  the  noise  fluctuation  process  and  (2)  some  of 
the  noise  dependency  is  accounted  for  in  selecting  p|.  On  the  other  hand,  the 
SNR  for  slow  fluctuations  is  more  sensitive  to  signal  power  fluctuations  than 
to  noise  power  fluctuations. 
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Appendix  A 
GA1MA  DISTRIBUTION 


The  normalized  Gamma  probability  density  function  (PDF)  is  the  standard 
Gamma  PDF  normalized  so  that  the  mean  is  independent  of  the  degrees  of  free¬ 
dom.  The  normalized  Gamma  PDF  is 


f(x)  = 


M-1  -MX/X0 
x  e _ 

(xo/M)Mr(M) 


x  >  0 


x  <  0 


(A.  1 ) 


where  N  is  the  degrees  of  freedom  and  x.  is  the  mean.  The  ath  moment  of  x  is 


M  =  E(x") 

a 


„  .  -Mx/x^ 

x*4*®"^  0  dx 


(xQ/M)"r(M)  Q 


HM±flIxa 
MaT(M)  ° 


(A. 2) 


Therefore,  the  mean  and  variance  are: 


x  =  Mi  =  xQ 


(A. 3a) 


M2  -  =  Xq/M 


(A. 3b) 


Note  that  the  variance  vanishes  as  M+«o. 


Let  x  and  y  be  two  independent  Gamma  distributed  random  variables  with 
degrees  of  freedom  Mg  and  My,  respectively,  and  mean  x0  and  yQ,  respec¬ 
tively.  Define  the  random  variable 


z  =  x/y 


(A. 4) 
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What  is  the  PDF  of  z?  Define  the  auxiliary  variable  w  =  y.  Then  the  Jacobian 
of  transformation  from  (x,y)  to  (w,z)  is 


J(x,y) 


i/y 

o 


Then, 


1 /y  =  1/w 


f(z,w)  =  wfxy(zw,w) 

and 


f(z)  = 


/<*>  < 

f(z,w)  dw  =  j 


wf  /xw,w)  dw 


Therefore, 

f(z) 


V1 


(V",f‘  r<Mx>  (wf 

/>v 


3  nwy) 


-(Mxl/*o  *  H  /y0)w 
e  dw 


M  -1 
x 


r<W  Mx 

l(M  )I  (M,  )  0  M  +M 


x  y' 


(1+az) 


x  y 


z  >  0 


where 


,  z  <  0 


a  = 


Vo 

Vo 


(A. 5) 


(A. 6a) 


(A. 6b) 


38 


v 

V 


The  mean  and  variance  of  z  are: 


x. 


2  "  V  yQ 


for  M  >1 

y 


4  =  (Mz>2  for  My  >  2 

x  y 


Define  the  random  variable 


z+1 


What  is  the  PDF  of  p?  It  is  easily  shown  that 


f(p)  = 


fz(1^ 

(1-p)2 


Substitue  Bq.  (A. 6a)  into  Eq.  (A. 9).  Then, 


f(p)  = 


r(M  +M  ) 

_ " _ jL _ 


r(Mx)  r(My) 


mx  A-WV1 


M  +M 

(Mo-Dp)  x  y 


,  0<P<1 


,  otherwise 


The  S'th  moment  of  p  is 


Mp(W 


E(pB) 

r(M  +m  )  ixm +8) 

f<wx)  r(Mx+My+By  2F1(MX+My’  Mx+B;  Mx+My+6;  1-<x) 


according  to  Eq.  3.197*3  of  Reference  A1 


«*.  **.  •  ■ '  V* 


(A. 7a) 


(A. 7b) 


(A. 8) 


(A. 9) 


(A. 10) 


(A. 11) 


Appendix  B 

EDGEWORTH  SERIES  FOR  COMPLEX  SPHERICALLY  BIVARIATE  RANDOM  PROCESSES 

The  derivation  of  the  Edgeworth  Series  for  complex,  spherically  invariant 
random  processes  is  presented  in  this  appendix.  A  random  process  is  spheri¬ 
cally  invariant  if  and  only  if  it  is  a  zero-mean  Gaussian  process  that  is 
multiplied  by  an  independent  random  variable  (Ref.  B.1).  The  derivation  of 
the  Edgeworth  Series  for  a  specific  type  of  spherically  invariant  random  pro¬ 
cess  is  obtained  by  (1)  computing  the  moment  generating  function,  (2)  comput¬ 
ing  the  cumulant  generating  function  from  the  moment  generating  function,  and 
(3)  finally  identifying  terms. 

B.1  Process  Description 

T 

Let  Z  =  (z1fZ2)  be  a  two-dimensional,  complex  zero  compound  process 
described  as 

where  S^  and  are  independent,  non-negative  random  variables  called  power 
processes,  which  are  also  independent  of  x^  and  y^;  xk  and  yk  are  independent, 
zero  mean  complex  Gaussian  random  variables  with  unit  variance;  ps  is  the  cor¬ 
relation  coefficient  of  x^  and  X2;  and  T  indicates  transpose.  Given  and 
N^,  Z  is  a  two-dimensional,  zero  mean,  complex  Gaussian  random  variable  with 
covariance  matrix 


where  *  indicates  complex  conjugate.  Therefore,  the  probability  density 
function  (PDF)  of  Z,  given  the  power  process,  is 


fCZ|s,,s2,Mt,M2>  = 


-z'irt 


(B.3) 


where  '  indicates  complex  conjugate  of  the  transpose.  Finally,  the  PDF  of  z 
is 


f(Z)  =  E(f(Z|SrS2,H,,N2)}S)S2il)iiN2 
where  u  jj  is  the  expectation  over  S-j  ,S2 ,N.j  ,N2< 


(B.4 ) 


B.2  Moment  Generating  Function 


The  moment  generating  function,  given  the  power  processes,  is 

<t>'R4> 


Mz(<(.|S1,S2,N1,N2)  =  e 


(B.5a) 


where 


♦  =  (4>-|  *<t>2^ 


(B.5b) 


Expand  Eq.  (B.5a): 


M2(4>|S1,S2,N1,N2)  =  exp  |  rn  |2  +  P2|<|>2|2  +  2v^S^  Re(ps4>*4>2)  } 


°°  i  &-k  k 

Sir  2  2Zcz,k,q,i 

Z=0  ksO  q=0  p=0 


.  i-k- 


*k-p+q.*  '  y*Z-k-q+p.-q+p 


(  R  .  ^ 


where 


r  -  fA-kuku  c  ^  2  nk_P  *<* 

C4,k,q,p  =  <•  q  Kp^lV  P1  p2  ps  ps 


(B.6b) 


and  (&)  is  the  binomial  coefficient.  From  the  discussion  in  section  B.1,  we 


know  that 


Mz(4>) 


=  E{Mz(*|S1,S2,N1,N2)}g^ ,S2,N1 ,N, 


,P4^) 


(B.7a) 


where 


A  A-k  k 


»»«♦>  ■ 

k=0  q=0  p=0 


(B.7b) 


^,k,q,p 


=  ^)(*'qk)tpK(k-P^-k)m2(p,l-k)pJ*k“qp*q  (B.7c) 


\(a.B)  =  E(rqS1®/2) 


B.3  Cumulant  Generating  Function 


(B.7d) 


The  cumulant  generating  function  is 


Kz(<|.)  =  An  (Mz(^)) 


=  An  (l  +  Mz ( 4> ) ) 


(B.8) 


where  M^O)  =  0.  In  light  of  the  discussion  in  section  (B.2),  Kz(0)  can  be 
expanded  in  two  ways.  First, 


2  VW 

r=1 


where 


(B.8a) 


E 

u+v+x+y=r 


-  uv*y  (»  )u(<t*)v(*  )x(<b*)y 
ulvlxly i  i  ^2; 


Ixly! 


(B.8b) 


is  a  polynomial  of  rth  order,  and 

gU+v+x+y 


uvxy 


3^  3^  3$. 


#y  z 

*2 


K  (♦) 


4b0 


(B.8c) 


The  second  way  is  to  substitute  Eq.  (B.7)  into  Eq.  (B.8): 

■  u 

.  .  i  I  n  t  M  ! 

co  (.!)' 


Kz'« 


,u»i  i  ye  fit) 


£cl 


u 


CB.9) 


Only  terms  up  to  r=4  in  Eq.  (B.8)  will  be  considered  because  of  the  com¬ 
plexity  of  the  problem.  By  expanding  Eqs.  (B.8)  and  (B.9)  and  identifying 
terms,  it  follows  that  Q's  are  related  to  the  P's  in  the  following  manner: 


en(  <*») 

=  03(<t>)  =  0 

(B. 10a) 

e2(  ♦) 

=  ?,(♦) 

(B.lOb) 

V# 

p*(<» 

•  P2<«  *  2 

e 

(B. 10c) 
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A>\  - . 


Also,  the  X«s  are  related  to  the  C*s  according  to  the  following 


e2($)\'a  : 

*2000  =  X0200  =  X0020  =  X0002 

X1010  =  X0101  3  0 

X1100  3  C1 100 

X0011  3  C1 101 

X0110  3  C1000 

X1001  3  C1010 

e4(*)X's  : 

X2200  3  4^2200  "  2Sl00 
—  —2 

X0022  3  4C2202  "  2C1 101 

—  —2 

X0220  3  4C2000  “  2C1000 

—  —2 

X2002  3  4C2020  “  2C1010 

X1210  =  2C2100  -  2C1000C1100 

X0121  3  2C2 101  "  2C1000C1101 

X2101  3  2C2 110  "  2C1010C1100 

X1012  3  2C2U1  "  2C1010C1101 


(B.  1 1 ) 


(B.12) 


(cont.) 


Therefore , 


Kz(^ 


e^*)  +  e2($) 

Siool^ll2  +  X001ll^2  +  W>2  +  X01104>14>2 


X2200i,  | 4  X0022ix  |4 

*  4  ^2 


*  V „ 

+  jj  V1  ?2  ♦  pf2 


►  2^r  1  V2  2  y vi  2 


*  XQ'3VA*  *  ^212U  Mf 

*  2*V2*2  2*1  2  2 


2ia  i  2 


B.4  Edgeworth  Series 


(B.14) 


Then, 


1100 

X01 10 

• 

0110 

X0011 

r11 

r12 

• 

r12 

r22 

.v*’ 

e 

(B. 15) 


(B. 16) 
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Substitute  Eq.  (B.14)  into  Eq.  CB- 16) : 


♦)  \®i,  (4*)J 

M  (*)  a  a  e  =  e 

z 


LJ~T\ 


i=0 

Using  only  the  first  two  terms  of  the  series  expansion  for  the 

e2(4>> 


Mz(<t»)  =  (1  +  94(4>))e 


.  X2200 i .  ,4  *0022..  ,4 

1  +  — — l<M  +-t“I*2I 


X0220.»2.2  X0220.2.»2 
+  ♦ 2  +  4  ^1^2 


X1210 .  .»2.  X1210 . 2 . 

♦  2  'Ml  *2  +  ~2  Vl*2 


X0121 . • .2. •  .  X0121 .  .  .»2 
-2-+1*2*2  +  ~2—V2*2 


It  is  easily  shown  that 


V>,R* 


_21_ 

2  e 

3rTl 


!♦  JV'W  = 


32  4'R4> 

2~  6 
3r22 


♦;%2e^’^ 


3  $'R$ 

—5-  eT 

3r12 


32  ♦,R<() 

.2  6 


23 


12 
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(B.17) 

exponential, 


(B.18) 

(B.19a) 

(B. 19b) 

(B. 19c) 

(B.19d) 


ViV* 


*11*12 


<fc,  *1  *2®  = 


*11*12 


<j>*  4>2$*e^' 
V1V2*2 


3r003r 


22  12 


^4*2  4*2  e 


*2_4>'R^  _  _ ai 


*22*12 


2U  , 2  ** R^>  _  *  6'K* 


|4»,  rl4>2re 


*1 1  *22 


Substitute  Eq.  (B. 19)  into  Eq.  (B. 18) -  Then, 


MU)  *  i  + 


*2200  32  \)022  S2 


4  Sr,2  4  *t2 


*0220  32  *0220 


4  *.2  4  .2 


*1210  32  X1210  32 

+  2  Sr.-Sr,,  *2  -  a  • 

11  12  *11*12 


*0121 


w 

\)121 


2  arTTorTT  2  ~  a  * 

22  12  3l223r12 


3r1 1 8r22 


(B.19e) 


(B.19f) 


(B.19g) 


( B .  1 9h ) 


(B.19i) 


(B.20) 


The  probability  density  function  of  z  is  obtained  by  taking  the  inverse 
Fourier  transform  of  Mz(0).  Taking  the  inverse  Fourier  transform  of 
Eq.  (B.20),  we  have 


g(Z) 


,  t  X2200  32  ^2200  32 

1  +  I.  «  +  .  . . . 


4  3n2 

3  11 

X0220  32 


4  3r2 

3r22 

\)220  32 


3r 


12 


3r 


12 


1210 


1210 


2  — 


3r  3r 
11  r12 


*0121  32 

+  ~T~W^ 


0121 


22  r12  2  3r  3r* 

22  r12 


+  A1 1 1 1  3r  3r 
11  r22 


f(Z) 


(f .21 ) 


where 


f(Z)  = 


1 


-Z'R-1Z 

e 


(B.22) 
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Appendix  C 

EDGEWORTH  SERIES  FOR  THE  CUMULATIVE  DERSITT  FUNCTION  OF 
THE  MSCC  FOR  A  SPHERICALLY  INVARIANT  PROCESS 

The  derivation  of  the  probability  and  cumulative  density  functions  (PDF 
and  CDF,  respectively)  of  the  sample  magnitude-squared  correlation  coefficient 
(MSCC)  is  difficult  for  non-Gaussian  signals.  The  Edgeworth  series  for  the 
PDF  of  the  sample  MSCC  will  be  developed  for  signals  that  are  spherically 
invariant.  The  approach  used  is  to  derive  the  PDF  of  the  sample  auto¬ 
covariance  matrix  of  the  observations  and  then  make  a  chnage  of  variables  to 
obtain  the  PDF  of  the  sample  MSCC.  The  processes  involved  in  the  derivation 
are  described  in  section  C.l.  The  derivation  of  the  characteristic  functions 
and  PDF  of  the  sample  auto-covariance  matrix  is  presented  in  sections  C.2  and 
C.3,  respectively.  The  PDF  of  the  sample  MSCC  is  obtained  from  the  PDF  of  the 
sample  auto-covariance  matrix  in  section  C.4.  Finally,  the  CDF  is  obtained  in 
section  C.5. 

C.l  Approach 

Let  Z&  be  a  two-dimensional  zero  mean  complex  random  column  vector  with 
elements  z-| (*•)  and  z2(M  representing  samples  from  channels  1  and  2  at  time 
*Tg  for  Z  =  1,2,...,Nt.  Ts  is  the  sampling  interval,  and  T  =  N^Tg  is  the 
observation  interval.  The  cross-covariance  matrix  of  Z(A)  is  defined  as: 

R  (l,k)  =  E(Z(4)  Z’(k)>  (C.l) 

z 

where  e(*}  denotes  statistical  expectation  and  '  is  the  complex  conjugate  of 
the  transpose.  Let 

zt  =  /s(A)  ♦  /nTTT  y4  (C.2) 

where  S(t)  and  N(l)  are  the  independent  two-dimensional  power  vectors  for 
signal  and  noise,  respectively;  X(X)  is  a  two-dimensional  unit-variance, 


zero-mean,  complex  Gaussian  random  vector  with  P^e  3  as  the  correlation 
coefficient  between  Xjd)  and  x2(4);  and  Y(Z)  is  a  two-dimensional  unit- 
variance,  zero-mean,  complex  Gaussian  random  vector  with  independent  compo¬ 
nents. 


Given  Sjj  and  N^,  is  a  two-dimensional,  zero  mean,  complex  Gaussian 
random  variable  with  covariance  matrix 


'  S1  +  N1 

0, 

R  = 

/S1S2  »! 

s2  +  N2 

r1 

sj:  p* 

1  2  Ks 

r2 

J 

where  •  indicates  complex  conjugate. 


The  exact  form  of  the  PDF  of  Zg  is  unknown.  However,  the  Edgeworth 
aeries  form  of  the  PDF  is  known  (see  Appendix  B).  The  PDF  of  is 


«(V 


(1  +  P)  f(Z^) 


(C.3) 


where 


P  s  — ! 


JL  C_2_£_ 

a_2  +  4  an2 

3rl1  3r22 


.i 

*  4 


2  3 


3r 


12 


.  9r„8r,2 


P! 


s  3r„JrJ2 


(cont.) 


2  La3r223r,12  S3r223rl2j 

C1  a  2^(2, 0)  -  B^I.O)2) 

C2  =  2(1112(2,0)  -  m2(1,0)2) 

C3  a  2(a1(0,2)  m2(0,2)  -  m^O.I)2  m2(0,1)2) 

C4  a  2 («1  (1,1)  b^O.O  -  m^(1 ,0)  Bj (0 , 1 )  82(0,1 ) ) 

Cg  a  2(B1(0,1)  Bgd.l)  -  Bgd.O)  B^O.1)  B^O.I)) 

C6  a  B^ (0,2)  B2(0,2)  -  B^o.l)2  B2(0,1)2 

Bk(a,&)  =  E{r®S^/2} 

1  "ZPZ£ 
f(Zft  S  n2 |r| 1/2  6 


E(r1  > 


PS 


E^S1S2>  PS 


E{r2) 


r  p 
11  12 


P12  r22 


The  aaaple  MSCC  can  be  coaputed  froa  the  sample  auto-cor relation 
aatrix.  The  two-dimensional  positive  definite  Hermetian  sample  auto¬ 
correlation  matrix  is 


where 


hZi 


Let 


A  = 


11 

L  a*2 


‘12 


22 


(C.10) 


(C.11) 


The  sample  MSCC  is  the  sample  magnitude-squared  cross-correlation  coefficient 
between  zj(&)  and  z2(l)  and  is  given  by 


I*, 2 1 


a  a 
1 1  22 


(C. 12) 


The  PDF  of  p2  can  be  derived  from  the  PDF  of  A  by  (1)  performing  the  change  of 
variables  indicated  in  Eq.  (C.6)  and  (2)  integrating  out  the  auxiliary  varia¬ 
bles  A] f  t  A22,  and  the  phase  angle  of  a12« 


C.2  Characteristic  Function  of  A 


The  characteristic  function  of  Aj,  using  the  Edgeworth  form  of  the  PDF  of 
Zt,  is 


%e>  . 

,  E  | 

s  11  *  P) 

li  -  jrJ-I 

nt 


(C.  13) 


where  |l  -  is  the  characteristic  function  of  A^  assuming  the  PDF  of 

is  f(Z^)  as  defined  in  Eqs.  (Ctf)  through  (C.8), 


Substitute  Eq.  (C.4)  into  Eq.  (C.13)  and  carry  out  the  indicated  partial  dif¬ 
ferentiation.  After  some  tedious  algebra, 


M.  (♦)  =  D— 1  ( 1  +  %) 

l  u 


where 


D  = 

Q  = 

a,  ■ 


li  -  j  *  fl 

T 

Q1  +  °2  +  C3 

\  b.  r2  +  i  p2 

I  2  22  2  r* 


11 


+  -^(p2r#2  ♦  p*2r2  ) 

+  2  l‘psr12  +  ps  r12J 

-  C„(p  r,_r„  ♦  p  r,_r__) 
4  Ks  12  22  s  12  22 


-  Wl2P11  +  pIri2r22) 


+  C6lps>  Criir22  +  ,ri2>  3 


ill 1 
sH 


°2  * 


-  cir22^y\  ~  C2rl1  *22 


♦  C^fp  r*  6*  4>  0*2r  A,  ^  ) 

3kMS  12V12  MS 


+  C4^Psr12*11  +  Psr12*11  "  tPs^12  r22 


Ps*12r22^ 


(C. 14a) 

(C. 14b) 
(C.l4c) 


(C.14d) 


(cont. ) 


*  C5(p.1*22  *  S2ri2*,1  -  ‘W„  *  P;*,2ru)) 

■  C6 ^ °s ^  ^11*11  *  r22*22  ‘  rl2**2  "  r*2*12^  {  "^3 

)  T 


(C.l4e) 


^  $2  +  ^  <(»2 
2  112  22 


+  C  (P2$«2  +  P»2<fr2  ) 
3  s  12  s  i2* 


+  C,  (P  4>  ♦  ■*■  p»4>  4>  ) 

4s  12  11  s  12  11 


♦  "  P3**,2*22> 


+  Cg  I  Pa  I  <*^*22  +  ^12^  *  1  „2 


(C.l4f) 


Since  the  Aj,'s  are  independent, 


M.(*)  *  |M.  (♦)]  T 


(C. 15) 


Substitute  Eq.  (C.l4a)  into  Eq.  (C.15): 


Ma(4>)  =  D  ^[l  ♦  Q/D2]1*1 

“^T  f  2 1 

~  D  |l  +  NTQ/D  | 


■  D  ♦ 


r(V2) 


(C. 16) 


C.3  Probability  Density  Function  of  A 

The  PDF  of  A  is  obtained  by  taking  the  inverse  Fourier  transform  of  MA( ♦) 
given  in  Eqs.  (C.14)  -  (C.16).  The  inverse  Fourier  transform  of  D~N  is  the 
complex  Uishart 


f(A|NT)  *  C(Nt)  V(a|nt) 


(C.17a) 
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where 


v(a|nt) 

C(Nt)  = 


M-  -N-TR(AS) 

1*1  • 


2NT,r,NT 


HJSJ _ , 

NT)r(NT  -  i 


nfoi 


’11 


3* 

12 


’12 


22 


(C. 17b) 


(C.  17c) 


(C.17d) 


which  was  derived  by  Goodman  (ref.  C.1).  Correspondingly,  the  PDF  represented 
by  the  characteristic  function  d~^+^  is  f(A|NT+2). 

-Nt 

By  making  use  of  the  fact  that  the  inverse  Fourier  transform  of  $akD 
can  be  obtained  by  (3<x/3a^lc)f(A|NT),  the  PDF  of  A  can  be  obtained  from  Mft($) 
by  taking  the  proper  partial  derivatives  with  respect  to  the  a?k  of 
f(A|Nf+2) .  Therefore,  substitute  Eqs.  (C. 17)  into  Eqs.  (C.ld)  -  (C.16)  and 
perform  the  indicated  partial  differentiation.  After  much  tedious  algebra, 


g(A)  =  fUlUj)  +  NtU(  A|  Nt)  (C.  18a) 

where  0  is  the  correction  term  of  f(A|NT)  which  is  the  PDF  of  A  for  complex 
Gaussian. 


U(A|Nt)  = 

U1(A|Nt) 

♦  U2(A|Nt) 

+  ^(AlHj) 

mJu1(a|nt) 

■  {> 

2  C2  2 
22  *2  11 

C3f  2  .2  .2  2  . 

♦  — (psr*2  * »;  ri2) 

-  <ps>-?2 

*  p;r,2)(C4 

r22  *  c5r11> 
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C.4  Probability  Density  Function  of  the  Sample  MSCC 

The  PDF  of  the  sample  MSCC,  p2,  is  obtained  from  the  PDF  of  A  by  (1)  per¬ 
forming  the  change  of  variables  indicated  in  Eq.  (C.12)  and  (2)  integrating 
out  the  auxiliary  variables  a^,  a22»  and  the  phase  angle  of  a^.  Performing 
this  operation  on  Eq.  (C.l8a),  we  have 

g(p2)  =  f(p2|NT)  +  NTf(p2|NT)  (C. 19) 

where 

_  _  N_  _  N_-2 

f(p2|NT)  =  (Nt  -  1 ) C 1  -  pp  (1  -  P  ) 

gF1 (Nt,  Nt;  1;  p2p2)  (C.19a) 

is  the  PDF  of  p2  for  no  fluctuations  and  Gaussian  signals  (Ref.  C.2);  f(p^|NT) 
is  the  correction  term  resulting  from  the  Edgeworth  series;  and 
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(e{^s^»2|p312 

(E{r1})2(E{r2»2 


(C.  19b) 


Since  f(P^lN^)  and  g(P^|N^)  are  PDF’s,  fCP2)^)  must  integrate  to  zero. 

Perform  the  indicated  change  of  variables  on  U(a|nt)  in  Eq.  (C.18)  to 
obtain  ftP2!^).  After  much  tedious  integral  evaluation  and  algebra, 
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C.5  Cumulative  Penalty  Function  of  the  Sample  MSCC 


The  CDF  of  the  sample  MSCC  is 


g(p^|nt)  =  /  fcG(  C*2  |ht)  d(^ 

o 

where  0  _<  p|  <:  1  is  the  threshold. 

Substitute  Eq.  (C.19)  into  Eq.  (C.23).  Then, 
g(^|ht)  =  f(P^|nt)  ♦  mtf(p^|nt) 


(C.22e) 


(C.22f) 


(C.22g) 


(C.23) 


(C.24) 


v'.  '.'..'Vw  '**  *\  I*j 


(C.25) 


F(P^|Nt)  =  (l-Pj)"1  L)(1-Pt)Vl(NT'k+1;1;PtPT) 

is  the  CDF  of  p2  for  no  fluctuation  (Ref.  C.2);  and 
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fPt 

FCp^JNt)  =  /  f(p*  |nt)  dt 
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Define 


f(p2|a,B,e,*,Y)  =  p2a(1-p2)62F1(0^;y;p2p2) 


(C.26a) 


and 


F(p2|a,B,e,*,Y) 


/ 

Jo 


f(p2|a,B ,6 ,$,y)  dp2 


(C.26b) 


F(p||NT)  has  the  same  form  as  f(p2|NT),  Eq.  (C.20),  where  the  function 
f(p2)a,B.0t$,Y)  is  replaced  by  F(p2|a,B,8,$ ,y) • 


Substitute  Eq.  (C.26a)  into  Eq.  (C.26b)  and  expand  the  hypergeometric 
function.  Then, 
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where 


I*(x+n) 

T(x) 


is  Pochanaier's  symbol.  From  reference  (C.3), 
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Therefore, 

?(Pt!NT)  =  PRiF(ptl°*MT“2,NTNT,1>  +  PR2F(Ptl0,NT'2,NT,NT+1,1) 

+  PR3F(p2|0,MT-2,HrNT+2,1)  ♦  PR4F(p2|0,NT-2,NT+1 ;NT+1 , 1 ) 

+PR^.F(p2|0,N^-1 ,Nt+1 ,Ht+1 ,1 )  ♦  PR6F(P2|0,Nt-1,Nt+1,Nt+2,1) 

+  pr7f(p2|o,nt,nt+2,mt+2,i)  ♦  pr8f(p2|i,nt-2,nt+i,nt+i,D 

♦  PRgF(P2|l,NT-2,NT4.1,NT+1,2) 

♦  PR10F(p2|l,NT-2,MT+2,NT+2,2) 

♦  PRnF(P2|l,NT-2,NT+2,NT+2,2) 

♦  PR12F(P2|2,Nt-2,Nt+2,Nt+2,3)  (C.30) 


where  the  PR's  are  defined  in  Eqs.  (C.21)  -  (C.22).  It  is  easily  shown  that 
?(1  |Hj)  «  0.  Therefore,  0(1 1%)  «  F(l|NT)  *  1,  which  makes  G(pf|NT)  a  CDF. 
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Appendix  D 

cosfpicibvt  evaluation  ns  thb  cop  of  the  sample  mscc  dnder 
rapid  fluctuation  CONDITIONS 

The  expression  for  the  Edgeworth  series  correction  factor,  F(Pt2(NT),  to 
the  cumulative  probability  function  (CDF)  of  the  sample  MSCC  contained  seven 
constants  (CN's)  which  depend  on  the  statistical  properties  of  the  fluctuation 
model,  Eqs.  (C.5)  and  (C.22)  of  Appendix  C.  These  constants  will  be  evaluated 
for  the  statistical  fluctuation  model  presented  in  Chapter  2.  The  important 


constants  are: 

C1  =  2(m1(2,0)  -  m,(1,0)2)  (D.la) 

C2  =  2^b2^2»°)  “  ^(I.O)2)  (D.lb) 

C3  =  2(m1(0,2)m2(0,2)  -  8^(0, 1)^(0, 1)2)  (D.le) 

cj|  =  2^(1, 1)^(0, 1)  -  m^l.Ojm^O.DmgtO.I))  (D.ld) 

C5  =  2(m1(0,1)m2(1,1)  -  m2(1,0)m1 (0,1 JfflgCO,!))  (D.le) 

C6  =  01^0,2)182(0,1)  -  m^O.D^O, I)2  (D.1f) 

where 

^(0,0)  s  E{r£sk0/2>  (D.2) 


rk  =  \  +  \  •  (D.3) 

Assume  that  is  Gamma  distributed  with  mean  and  MS^  degrees  of  free¬ 
dom  and  N|(  is  also  Gamma  distributed  with  mean  %  and  MN^  degrees  of  freedom 
as  disoussed  in  Chapter  2.  Then,  from  Eq.  (A. 2)  of  Appendix  A, 
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According  to  Eqs.  (2.17)  and  (2.18)  of  Chapter  2, 
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where 
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(D.6b) 


is  the  bias  factor,  and 

SNRk  = 

vNk  . 

(D.6c) 

Substitute  Eqs.  (D.5)  and  (D.6)  into  Eq.  (C.20)  of  Appendix  C. 

The  CN’s 
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